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We present a comprehensive inelastic neutron scattering study of the magnetic excitations in twin-free 
YBa 2 Cu 3 6 6 (T c = 61 K) for 5 K < T < 290 K. Taking full account of the instrumental resolution, we derive 
analytical model functions for the magnetic susceptibility /"(Q, co) at T = 5 K and 70 K in absolute units. 
Our models are supported by previous results on similar samples and are valid at least up to excitation 
energies of co = 100 meV The detailed knowledge of /"(Q, co) permits quantitative comparison to the 
results of complementary techniques including angle-resolved photoemission spectroscopy (ARPES), as 
demonstrated in Dahm et ah, Nature Phys. 5, 217, (2009). Based on accurate modeling of the effect of 
the resolution function on the detected intensity, we determine important intrinsic features of the spin 
excitation spectrum, with a focus on the differences above and below T c . In particular, at T = 70 K the 
spectrum exhibits a pronounced twofold in-plane anisotropy at low energies, which evolves towards fourfold 
rotational symmetry at high energies, and the related dispersion is "Y"-shaped. At T = 5 K, on the other 
hand, the spectrum develops a continuous, downward-dispersing "resonant" mode with weaker in-plane 
anisotropy. We understand this topology change as arising from the competition between superconductivity 
and the same electronic liquid-crystal state as observed in YBa 2 Cu 3 6 4 5 . We discuss our data in the context 
of different theoretical scenarios suggested to explain this state. 
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I. INTRODUCTION 



In the recent years, substantial progress has been made in 
the characterization of fermionic quasiparticles in metallic 
and superconducting cuprates. While quantum oscillation 
experiments on both overdoped and underdoped cuprates 
document the presence of Landau quasiparticles on at least 
some segments on the Fermi surface, 1-4 a variety of spectro- 
scopic experiments including angle-resolved photoemission 
spectroscopy (ARPES), 5-8 scanning tunneling spectroscopy 
(STS) 9 ' 10 , and optical spectroscopy 11 ' 12 demonstrate that 
these quasiparticles are strongly renormalized by interac- 
tions with low-energy bosonic modes. A detailed descrip- 
tion of this fermion-boson coupling may hold the key to 
a quantitative understanding of high-temperature super- 
conductivity. However, in order to assess the viability of 
boson-exchange pairing mechanisms in general 13 ' 14 as well 
as the strength of particular electron-boson pairing chan- 
nels, the bosonic spectral functions must be determined as 
accurately as the fermionic ones. For spin excitations and 
phonons, the bosonic excitations currently under discussion 
as the main candidates for strong coupling to fermions, this 
can be accomplished by inelastic neutron scattering. The 
neutron scattering data can then be cross-correlated with 
the results of complementary probes that provide direct 
information on the fermionic excitations. 



Here we report a comprehensive inelastic magnetic neu- 
tron scattering study of spin excitations in the underdoped 
high-temperature superconductor YBa 2 Cu 3 6 6 (hole con- 
centration p h ~ 0.12, superconducting transition tem- 
perature T c = 61 K), which is suitable as a model com- 
pound for several reasons. First, unlike the widely stud- 
ied Bi 2 Sr 2 Ca n _ 1 Cu n 02( n+ 2)+5 system, the crystal structure 
of YBa 2 Cu 3 6 6 does not include incommensurate lattice 
modulations. Moreover, its electron system is less af- 
fected by chemical disorder than those of many other high- 
temperature superconductors, as demonstrated by excep- 
tionally narrow nuclear-magnetic-resonance lines, 15 and 
by the fact that quantum oscillations were observed in 
YBa 2 Cu 3 6+JC crystals with somewhat lower doping levels. 1 
Crystallographic twin boundaries, a common type of mi- 
crostructural defect that disturbs the electron system and 
obscures detection of anisotropics in the Cu0 2 planes, can 
be effectively eliminated in YBa 2 Cu 3 6 6 . Finally, the doping 
level of YBa 2 Cu 3 6 6 is high enough such that the effect of 
superconductivity on the spin dynamics is clearly apparent, 
but low enough to ascertain that the spin excitations are suf- 
ficiently intense to allow neutron scattering measurements 
with high signal-to-background ratio. 

A series of prior neutron scattering studies has already 
provided valuable information about the overall layout of 
the spin excitation spectrum of YBa 2 Cu 3 66 . 16-18 In par- 
ticular, these studies showed that the spectrum in the su- 
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perconducting state exhibits the so-called hourglass disper- 
sion in which upward- and downward-dispersing branches 
meet and form the commensurate resonance peak at the 
wave vector characterizing antiferromagnetism in the un- 
doped parent compounds. Similar hourglass dispersions 
were also reported for optimally-doped YBa 2 Cu 3 6+x , 19 ' 20 
and other families of hole-doped cuprates including 
La 2 _ x (Sr,Ba) x Cu0 4 , 21 ' 22 and Bi 2 Sr 2 CaCu 2 8+5 . 23 ' 24 How- 
ever, since the small neutron scattering cross section requires 
crystals with volumes in the cm 3 range, most of the earlier 
experiments on YBa 2 Cu 3 6 6 were carried out on specimens 
that contain parasitic phases as well as microcracks and 
other microstructural defects, and are hence not suitable 
for measurements with most other experimental methods. 25 
Since samples grown in different laboratories often vary in 
their physical properties even if their compositions are nomi- 
nally identical, this situation has severely limited any quanti- 
tative comparison of neutron scattering data to the results of 
complementary probes. In order to overcome these difficul- 
ties, we have performed neutron scattering experiments on 
mosaics consisting of small YBa 2 Cu 3 6 6 crystals of superior 
quality. Specimens grown under identical conditions have 
turned out to be suitable for investigation by a variety of 
complementary experimental techniques 5-7 ' 18 ' 19 ' 26-35 . 

Following this strategy, we have recently taken the 
first steps towards a quantitative cross-correlation of neu- 
tron scattering and ARPES spectra measured on the same 
YBa 2 Cu 3 6 6 crystals. 33 An analytical fitting function for the 
imaginary part of the wave-vector and energy dependent 
spin susceptibility, #"(Q, co), determined by neutron scat- 
tering proved to be an essential ingredient in this study. 
We expect this approach to be more generally useful for at- 
tempts to determine the coupling between spin excitations 
and fermionic quasiparticles from spectroscopic data, and 
to asses its role in the formation of the superconducting 
state. In this article, we hence provide a comprehensive 
description of the spin dynamics of YBa 2 Cu 3 6 6 at temper- 
atures above and below T c . We provide a full account of 
the influence of the instrumental resolution function on the 
intensity and shape of the features observed in the neutron 
spectra, and we test various fitting models against the full 
experimental data set. The outcome is an analytical function 
for x"(Q, co) in absolute units that reliably reproduces all 
available experimental data up to at least co ~ 100 meV 

The second major motivation for our study derives from 
theoretical work that has been inspired by the hourglass 
dispersion. The scenarios proposed to describe this unusual 
excitation mode range from itinerant-electron models in 
which the hourglass is due to an excitonic collective mode 
of a d-wave superconductor, 36 ' 37 to collective excitations of 
a variety of ordered states with localized electrons that com- 
pete with d-wave superconductivity, including site-centered 
and bond-centered stripes, 38-42 , spin spirals 43-46 and models 
which take into account both itinerant excitations and lo- 
cal spin rotations. 47 These competing states break both the 
translational and the rotational C 4 symmetry of the Cu0 2 
square lattice and thus generate Bragg reflections in neutron 
scattering experiments. Bragg peaks characteristic of incom- 



mensurate magnetic order have indeed been observed at 
low temperatures in La 2 _ x (Sr,Ba) x Cu0 4 in wide range cen- 
tered around x ~ p ~ 1/8, 48-51 and in deeply underdoped 
YBa 2 Cu 3 6+x . 31 The fact that no such peaks are present in 
YBa 2 Cu 3 6 6 indicates that for this composition fluctuating 
stripe 52 or fully itinerant approaches 36 ' 53 might be a better 
starting point. 

Recent neutron scattering experiments on untwinned 
YBa 2 Cu 3 6+x crystals with a somewhat lower doping level 
(x = 0.45, p h ~ 0.085) have provided evidence of a col- 
lective ordering phenomenon besides static magnetic order 
and d-wave superconductivity. At this doping level, static 
magnetic order vanishes upon heating above ~ 2 K, but 
low-energy incommensurate spin fluctuation are observed 
up to much higher temperatures. Since the incommensura- 
bility selects a unique axis in the Cu0 2 plane and exhibits a 
sharp onset at T ~ 150 K, these data have been interpreted 
in terms of an "electronic liquid crystal" phase that breaks 
the C 4 rotational symmetry down to C 2 , but conserves the 
translational symmetry of the lattice. 54 This phase can be 
viewed from a localized-electron perspective as an array of 
fluctuating stripes or spirals with a spontaneously selected 
propagation vector in the Cu0 2 planes, 52 ' 54 ' 55 or from an 
itinerant-electron perspective as a Pomeranchuk state with 
a spontaneous Fermi-surface deformation. 56-60 The small 
orthorhombic distortion of the YBa 2 Cu 3 6+x crystal struc- 
ture serves as an aligning field for the propagation vector 
that permits direct measurements of the in-plane anisotropy 
with volume-integrating experimental probes. Strongly tem- 
perature dependent in-plane anisotropics of the electrical 
conductivity 61 ' 62 and the Nernst effect 63 are also consis- 
tent with electronic liquid-crystal behavior in YBa 2 Cu 3 6+x , 
as are neutron scattering results for strongly underdoped 
L^i .96 Sr 04 CuO 4 . 64 

In a recent brief communication, 18 we had shown that 
the spin dynamics of YBa 2 Cu 3 6 6 at temperatures above 
T c also exhibits a strong in-plane anisotropy reminiscent of 
electronic liquid-crystal behavior, consistent with the obser- 
vation of strongly anisotropic charge transport at this doping 
level, and that the crossover to the more isotropic dynam- 
ics below T c involves a change in the topology of the spin 
excitation spectrum. Here we fully document this topology 
change, with particular focus on the effect of the instrumen- 
tal resolution function. We expect that the comprehensive 
description of the spin dynamics above and below T c , in- 
cluding the analytical fitting functions, will prove useful as 
a basis for comparison to theoretical work on the spin dy- 
namics of electronic liquid crystals 45 ' 46 ' 52 ' 65-70 and d-wave 
superconductors 36 ' 71-78 , and on the interplay between these 
two ordering phenomena. 

The paper is structured as follows. In section II we de- 
scribe the experimental method and technical details. In 
II A we define the physical quantities we have measured and 
the notation we adopt. In II B we discuss the experimental 
setup, with particular emphasis on aspects relevant for the 
evaluation of the instrumental resolution. We further ex- 
plain the choice of the scattering geometry. In II C we give 
details of the data analysis, the background correction, and 
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the way different parts of the data measured under different 
conditions were brought to the same scale. The characteris- 
tic properties of the YBa 2 Cu 3 6 6 sample we have used are 
summarized in II D. In section III we show our experimental 
results. First, in III A we give a qualitative overview of the 
most salient features of the spectrum, already apparent in 
the raw data. Then we present data the below (III B) and 
above (III C)) the resonance energy co T ~ 38.5 meV in more 
detail, both for the superconducting and the normal states. 
In III D we discuss the dispersion relations in both states. In 
III E, finally we show the %"(Q, co) models which we found 
to fit best our data in the superconducting and the normal 
state, respectively. At the end of III E we discuss the validity 
range of our models and compare to the data previously 
reported. In section IV, we discuss the implications of our 
data for different proposed states with broken C 4 symme- 
try and the connection to the electronic liquid-crystal state 
observed in more strongly underdoped YBa 2 Cu 3 6 45 . 



II. EXPERIMENTAL METHOD AND TECHNICAL DETAILS 

A. Neutron scattering definitions and notation 

The differential neutron scattering cross section is related 
to the spin-spin correlation function S(Q, co) in the following 
way (in the general case of magnetic contributions to the 
scattering, S(Q, co) is a tensor): 79 



s 

af3 



5 a/3 ~ 



S ap (Q,co). (1) 



Here, r is the classical electron radius, y = 1.913, N is 
the number of spins in the system, k { and k f are the wave 
vectors of the incident and scattered neutron, Q = k f — k { , 
and a and f5 are the components of S and Q. Further, we 
adopt the notation that co is the energy lost by the neutron 
(in our notation, ft is 1). The Debye-Waller factor e ~ 2W ^ 
tends to suppress the intensity at higher Q, but within the 
error it is not important at the temperatures we work and 
will be neglected. The magnetic form factor /(Q), which 
represents the Fourier transform of the unpaired electrons 
in orbitals contributing to the scattering (in our case, in the 
Cu 2+ x 2 — y 2 orbital), suppresses the intensity at higher Q 
as well. It is taken into account as described in II C. 

The imaginary part of the magnetic susceptibility tensor 
Xap(Q> to) 5 which defines the response of the system, is 
the more straight-forward quantity for theoretical consider- 
ations. Hence, the model functions derived in III E represent 
this quantity. Susceptibility and correlation function are 
related by 

x'a P (Q,<») = <£M B ) 2 [1 " exp(-tol3)]S ap (Q,co) (2) 

where f5 = l/k B T. At the temperatures for which we de- 
velop our analytical models in III E, the detailed balance 



factor 1 + n(co) = 1/[1 - exp(-co/3)] can be neglected 
above ~ lOmeY i. e. essentially for all energies we have 
studied, and therefore, S oc %" . 

In the absence of magnetic long range order and a pre- 
ferred spin orientation, Y, af5 {5 a§ - Q a Qp/Q 2 )Xap^ to) = 
2Tr (*ais)/ 3 - Usin S the isotro P ic susceptibility %" = Tr(^), 
the scattering cross section becomes 79 ' 80 
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Here, we adopt the notation that co is the energy lost 
by the neutron (we set ft to 1), while Q is the wave- 
vector transfer to the neutron. Q = (H,K,L) is given in 
units of the reciprocal lattice vectors a*, b* and c* (r.l.u., 
as usual omitted in the following), where a = 2n/a* = 

3.82A" 1 , b = 2n/b* = 3.875A" 1 and c = 2ti/c* = 

o -1 

11.73 A . As in the parent compound YBa 2 Cu 3 6 , mag- 
netic excitations are centered around the antiferromagnetic 
wavevector Qaf M = (0.5, 0.5, L), and the equivalent points 
([2ft + 1] 0.5, [2k + 1] 0.5, L) in higher Brillouin zones. 

The YBa 2 Cu 3 6+JC structure comprises two antiferromag- 
netically coupled Cu0 2 layers separated by d = 3.285 A 
(bilayers), while the distance between adjacent bilayers is 
much larger and the coupling is three orders of magnitude 
smaller. Using a Heisenberg model for the parent compound, 
it can be shown that the intra-bilayer coupling gives rise to 
acoustic and optical magnetic excitations, whose intensity is 
modulated by G 2 dd and G 2 ven , where G odd = sm^nLd/c) and 
G even = cos(nLd/c) are the acoustic and optical structure 
factors, respectively: 81 
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This bilayer modulation persists well into the overdoped 
regime, 19 ' 28 and even excitations exhibit a gap of about 
30 — 40meV in our doping range. 80 In this study we will 
concentrate exclusively on odd excitations, which peak at 
L ~ {21 + 1) x 1.7. The even channel has been previously 
studied in some detail, 19 ' 28 ' 80 ' 82 however the determination 
of an analytical formula describing #^ en (Q, co) needs to 
await future investigations. 



B. Inelastic neutron scattering setup 

The technique of three-axes neutron scattering is de- 
scribed in Ref. 83. The signal measured in a scattering 
experiment is proportional to the differential scattering cross 
section convoluted with the instrumental resolution func- 
tion. The resolution function depends in a complicated way 
on the properties of the instrumental setup. 83 In the fol- 
lowing, we will describe these properties and refer to II C 
for details about the software used to perform the actual 
resolution analysis. 

Most of the results presented here were obtained at 
the thermal three-axes spectrometer IN8 (Institute Laue- 
Langevin, Grenoble), while some measurements were done 
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Fig. 1 (color online). Experimental setup and sample characteristics, a) The "W-configuration" of the IN8 spectrometer at ILL, top view. 
The angle names follow the ILL nomenclature. The table shows the sizes of the components and distances between them and further setup 
details. All measures are in cm. The angles A2, A4 and A6 are defined to be zero if the incoming and outgoing beam at the corresponding 
component are parallel. The virtual source is defined by an "empty collimator" (w = 3 cm, h — 16 cm and I — 18 cm), placed between the 
real source and the monochromator. The distance between the latter two is 810 cm. b) Longitudinal scans through the (200) and (020) 
structural Bragg reflections demonstrating a twin-domain population ratio of 94:6. 



at the thermal spectrometer 2T (Laboratoire Leon Brilloiun, 
Gif-sur-Yvette, France) and at the time-of-flight (TOF) spec- 
trometer MAPS (ISIS spallation source, Rutherford Appleton 
Laboratory, Chilton, UK). 

In Fig. la) we show a sketch of the IN8 setup, along 
with characteristics of the components and the distances 
between them. Both the monochromator and the analyzer 
consisted of an array of pyrolytic graphite (PG) crystals 
set for the (002) reflection, which allowed both horizontal 
and vertical focussing to maximize the neutron flux at the 
sample. We worked in the constant-/c f mode, and a PG filter 
was put in the path of the scattered beam to extinguish the 
contamination from higher-order neutrons. This filter is 

particularly well suited for measurements at fc f = 2.66 A 

and is slightly less efficient at k { = 4.1 and 4.5 A . Further 
constraints to be taken into account are the kinematics of 
the neutron scattering process which restricts measurements 
at low fc f to low energy transfers, and the fact that the 
resolution is better at low fc f , usually at the expense of 
a lower neutron intensity. Consequently, measurements 

were done at fc f = 2. 66 A up to the resonance energy of 
co r = 38 meV which roughly coincides with the kinematically 
allowed range (in the given BZ, see below) and also the 
range of most intense excitations (Fig 4). Above co T , the 

highest reasonable k { of 4.5 A was used, again dictated by 
kinematic restrictions at high energies. 84 

A special feature of IN8 is its virtual-source configuration, 
employed to maximize neutron flux: The real source is circu- 
lar, has a large diameter of 20 cm and is 810 cm away from 
the monochromator. In order to use the horizontal focussing 
mode of the monochromator without degrading the energy 
resolution, the distances L and L x must be identical (Row- 
land condition), which is achieved by a slit appropriately 
placed between real source and monochromator, thus defin- 



ing the virtual source. In principle, the Rowland condition 
must also be fulfilled for the horizontally focussing analyzer, 
however one can choose L 3 smaller if the horizontal extent 
of the detector is large enough to accept all the incoming 
neutrons of the beam focussed behind the detector plane. 

A further requirement to fulfill the Rowland condition is 
the appropriate choice of the k { (fc f ) dependent focussing 
curvature to ensure that the Bragg angle is constant across 
the monochromator (analyzer) surface. Usually, the cur- 
vature is optimized empirically to obtain the best energy 
resolution at the instrument. For the resolution calculations 
(see also III E) we will use the following formulae: 

1 sinf3 M 1 sin0 A 

^MH ^1 ^AH ^2 

where R MH and R m are the horizontal curvature radii and M 
and A are the Bragg angles at the monochromator (MH) 
and analyzer (AH), respectively. 

Since to a first approximation the vertical focussing has 
no impact on the Bragg angle, only the vertical resolution is 
affected which, as usual, was quite relaxed and optimized 
for highest flux in our case. For the calculations we use the 
following formula for the vertical curvature radius: 



R v V^bef LaftJ 2 sin0 

where L bef = L or L 2 , L aft = L 1 or L 3 and = 9 M or 
A for the monochromator and analyzer, respectively. The 
vertical analyzer focussing is fixed and optimized for an 

intermediate fc f of ~ 3.3 A which was taken into account 
in the resolution calculations. 

Since we are interested in the spin excitations in the Cu0 2 
planes and hence need a high resolution there, it seems opti- 
mal to work in a configuration with the a*- and b* -directions 
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are in the horizontal plane. However, since we are probing 
excitations in the odd channel, whose intensity disappears 
at L = 0, and its first maximum appears for L ~ 1.7, see II A. 
Taking into account the further constraints - to keep Q small 
to gain intensity due to the Cu 2+ magnetic form factor and 
to put c* as close to vertical as possible - we measured most 
scans along a* and b* around Q AFM = (0.5, —1.5, —1.7) and 
(1.5,0.5, 1.7), respectively. This entails a scattering plane 
in which the c*-axis is tilted by ~ 19° from the vertical 
direction. 

An important objective of this study is to see how the 
excitations evolve around Q AFM in the two distinct directions 
a* and b*. Thus, to eliminate any remaining uncertainties 
related to details of the resolution calculations, we decided 
to perform the scans in these two directions under identical 
resolution conditions. This was achieved by performing the 
a*- and b* -scans around the two equivalent Q AFM indicated 
above. This change of scattering planes corresponds to a 
rotation of the resolution ellipsoid by 90° w. r. t. the a-b- 
plane. To avoid the necessity of taking the sample out of 
the cryostat for this change, we developed a special holder 
with two cradles which allowed us to "pre-tilt" the sample 
by ~ 10° around the two perpendicular directions a* and 
b*, respectively, before mounting it into a standard orange 
cryostat. The cryostat itself, fixed to the sample table, can 
only be tilted by about ±10° around the same directions 
- our special holder effectively shifts the zero by 10° and 
the accessible tilting range to ~ 0° — 20° for both directions. 
This allowed us to access the necessary angle of 19° and 
change between the a*- and b* -configuration in a simple 
way. 



C. Data analysis 

At our doping level, the spectral weight of the spin excita- 
tions and the signal-to-background ratio are relatively large, 
so that at most energies a subtraction of a linear background 
provides sufficient correction (raw data are shown in Fig. 2) . 
At energies ~ 30 meV there are sometimes spurious peaks, 
which, however, turned out to be temperature-independent 
in most cases. Since at these energies the spin excitations are 
strongly overdamped close to room temperature (see also 
Fig. 2), the magnetic contribution to the neutron-scattering 
cross section is best represented by the data from which a 
background scan obtained at T = 250 K has been subtracted. 
For the sake of consistency, in several cases we perform this 
correction up to co T = 38 meV, e. g. in the color maps of Figs. 
8 and 9, or in the fits of Fig. 17. Naturally, due to the smooth 
background there (Fig. 2e-h) this has no consequences for 
our analysis. 

As mentioned, the energy ranges below and above co r 
were measured with different fc f and therefore under dif- 
ferent resolution conditions. Furthermore, kinematical con- 
straints forced us to measure scans at co > 70meV in 
higher Brillouin zones (Fig. 7) where the magnetic inten- 
sity is suppressed due to the Cu 2+ form factor (which is 
anisotropic 85-87 and only known with limited precision). 



For a quantitative analysis, it is thus important to put mea- 
surements done under such different conditions onto the 
same intensity scale. To this end, we measured the reso- 
nance peak at co r = 38 meV with all fc f and in all Brillouin 
zones used in these experiments. We then fitted the scan 
profiles to the damped harmonic oscillator (DHO) model 
of the hourglass dispersion discussed in HIE, taking full 
account of the instrumental resolution, and extracted corre- 
sponding scaling factors. This allowed us to obtain a very 
satisfactory fit of the measured scan profiles at all energies 
to the same DHO -model with the same, globally valid pa- 
rameter values (Figs. 17 and 18). This, and the fact that the 
energy evolution of the Q-integrated susceptibility of the 
model agrees well with the evolution of the Q-integrated 
susceptibility obtained directly from the data (Fig. 12) lends 
further support to the validity of our approach. 

We used the program Resca.15 88 to calculate the resolu- 
tion function and convolute it with the model functions 
we develop in III E. Our resolution calculations were based 
on the Popovici method 89 which takes full account of ge- 
ometrical properties of the instrument components (sizes 
and distances) . Since the program had not been used for 
setups with a virtual source before, we compared the ob- 
tained resolution function with the results of the program 
RESTRAX5.0.5 90 which uses ray-tracing methods and is thus 
more flexible in implementing unusual geometries. The dis- 
crepancy between the two programs concerning the resolu- 
tion width was below 10% and justifies the usage of Rescal5, 
which allows the import of data, definition of parameters 
and other things in a more transparent way. 



D. Characteristics of the YBa 2 Cu 3 6 6 sample 

The sample on which the experiments were carried out 
is an array consisting of 180 co-aligned YBa 2 Cu 3 6 6 single 
crystals, with a total volume of ~ 450 mm 3 and a total 
mosaicity of less than 1.2°. The crystals were detwinned 
by subjecting them to a uniaxial mechanical stress of ~ 
5 x 10 7 Nm" 2 along a (100) direction at 400 °C in argon 
for at least two hours. 91 ' 92 After detwinning, they exhibit 
superconducting transition temperatures (midpoint) of T c ~ 
61 K and transition widths AT C ~ 2K, as determined by 
magnetometry for each crystal. A comparison of the area of 
the main (200) and (020) Bragg reflections to the area of 
the respective satellites due to the minority twin domains 
yields a twin-domain population ratio of 94:6 for the entire 
sample, Fig. lb). The doping level is ~ 0.12 holes per planar 
copper. 93 ' 94 Further preparation details are given in Refs. 18 
and 26. 



III. RESULTS 

This section is structured as follows: In III A we give a 
brief qualitative overview of the most salient properties of 
the spectrum, already apparent in the raw three-axes (Fig. 2) 
and time-of-flight (TOF) data (Fig. 3), and the energy evo- 
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Fig. 2 (color online) . Energy evolution of the in-plane spin exci- 
tations around Qaf M at different temperatures. The raw data are 
shown at 5 K and 70 K, while a constant BG has been subtracted at 
250 K and 290 K to take into account the increase of multi-phonon 
scattering. Corrections for the Bose factor are small and were 
not applied. The final wave-vector k f was fixed to 2.66 A for 

co < 37meV and to 4.5A at higher energies. Scans along a* and 
along b* were performed under identical resolution conditions. 
The data were previously shown in Ref. 18. 



A. Overview 

In Fig. 2 we show unprocessed ("raw") scan profiles for 
different co throughout the covered energy range, measured 
along the two principal directions a* and b*. Several charac- 
teristic aspects of the data which we will investigate in quan- 
titative detail later on are visible here on a qualitative level: 
Beginning deep in the superconducting state at T = 5 K, we 
observe, as previously reported, that the incommensurate 
(IC) branches at high 16 ' 17 ' 95 and low 96 ' 97 energies merge to 
a single, commensurate resonance peak at co r ~ 38 meV The 
resonance peak, 98 originally discovered in YBa 2 Cu 3 6+x , 99 
was later also found in bilayer Bi 2 Sr 2 CaCu20 8+ 5, 100 single- 
layer Tl 2 Ba 2 Cu0 6+5 (Ref. 101) and HgBa 2 Cu0 4+5 , 102 and 
in electron-doped cuprates. 103 ' 104 The low- and high-energy 
branches, together with the resonance peak, form the "hour- 
glass" dispersion with commensurate neck. 16-18 ' 22 ' 95 Not 
only the peak intensity but also the Q-integrated intensity 
are maximal at the resonance energy co r (Fig. 4). 99 It is note- 
worthy that although there is a noticeable a- b -anisotropy 
at low energy, the signal is two-dimensional and incom- 
mensurate in both directions. Further properties become 
apparent from the constant-energy TOF maps in Fig. 3: First, 
the anisotropy increases with decreasing energy; second, 
the signal is nearly four-fold symmetric at high energies 
above co r ; 18,105 third, the orientation of the intensity max- 
ima at high energies is rotated by 45° with respect to the 
distribution at low energies. 17 ' 22 

Upon raising T to 70 K, i. e. 10 K above T c , prominent 



30 meV 



lution of the Q-integrated intensity (Fig. 4) . In III B and 
III C we discuss in more detail the in-plane anisotropy and 
temperature evolution at low energies (< co r , Figs. 5 and 
6) and high energies (> co r , Fig. 7), respectively. In HID, 
we first describe a novel way to visualize the dispersion of 
data with strongly energy-dependent intensity in energy- 
momentum maps (Figs. 8 and 9). We then compare the 
dispersion topology below and above T c both using these 
maps and a more conventional representation, namely plots 
of the energy evolution of the incommensurabilty A IC , as 
obtained from Gaussian fits to the scan profiles (Fig. 10 and 
19). Based on this comparison and on a careful evaluation 
of the data in the vicinity of co r , we show that the hourglass 
dispersion with commensurate "neck" observed in the su- 
perconducting state is replaced by a Y-shaped dispersion: 
The latter is very steep and quasi one-dimensional below 
~ co T (Fig. 11), and exhibits an incommensurate geometry 
at any energy, including co r . Finally, in III E we develop em- 
pirical, analytical model functions for #"(Q, co) based on a 
DHO model at T = 5 K, III E 1, and a combination of a DHO 
model and Gaussian profiles at T = 70 K, HIE 2. Despite the 
small number of free parameters, the models exhibit all the 
spectral features described in III A - III D and provide a very 
good fit to the experimental data (Figs. 12-18), with the 
same parameter values valid at all energies. 




0.2 0.4 0.6 0.8 
H (r.l.u.) 



Fig. 3 (color online). Color representation of the in-plane magnetic 
intensity collected at the TOF spectrometer MAPS at ISIS for dif- 
ferent energy transfers. The temperature was 3 K. Intensity was 
integrated over an energy range of ±3.5 meV in a) and b), ±5 meV 
in c) and ±7.5 meV in d). The data represent acoustic excitations 
and were collected at L = 1.6 in a), at 2.1 in b) and at 4.7 in c) 
and d), i. e. close to one of the maxima at L = 1.7 and 5.1. The 
incident energy was fixed to 120 meV in a) - c) and to 160 meV in 
d). The source proton current was 170 //A and the Fermi-chopper 
rotation frequency was set to 250 Hz. A background quadratic in 
Q was subtracted, and data are shown in arbitrary units. 
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Fig. 4 (color online). Imaginary part of the Q-integrated (local) 
susceptibility as a function of the excitation energy at 5 K, 70 K and 
290 K. 



changes occur at and below the resonance energy while 
the high energies are hardly affected. This supports the no- 
tion that the resonance peak and the downward dispersing 
branch form a common, superconductivity-induced reso- 
nance mode. 95 Notably above T c this mode is replaced by 
an incommensurate, nearly one-dimensional distribution. 18 
In the vicinity of co r and down to ~ 15 meV below co r , the 
Q-integrated intensity is higher in the superconducting state 
(Fig. 4), reflecting the appearance of the resonance mode. 
Fig. 4 also indicates the opening of a superconductivity- 
induced gap below ~ 20meV 106 Neglecting the redistri- 
bution in Q at constant energy for the moment, the reso- 
nance mode gathers its spectral weight mostly from this 
gap and from a limited energy range around ~ 45 meV 
Altogether, the impact of superconductivity is limited to 
co < 50 — 55 meV, and, concerning the excitation topology 
alone, even to co < co r (see also Fig. 8e,f). 

Finally, raising T to room temperature eradicates any 
sign of coherent excitations for co < co T and also notably 
suppresses the intensity at higher energies. At the same 
time, the incommensurability A IC decreases. 



B. Low-energy excitations 

In Fig. 5 we show Q-scan profiles in an extended energy 
range co < co r , below (5 K) and above (70 K) T c . The data 
in the superconducting state exhibit a continuous, mono- 
tonic dispersion which is somewhat steeper along b*, and 
exhibits a marked in-plane intensity anisotropy, 18 ' 26 as first 
reported in Ref. 107. When heating to 70 K, we observe 
marked changes in the entire energy range: First, the in- 
plane distribution of the intensity changes, and we observe 
a quasi-one-dimensional signal in all measurements. The 
data can be well fitted with two Gaussian peaks displaced 
in the a*-direction from Q A fm> whereas a single Gaussian is 
sufficient along b* and the fit is not significantly improved 
assuming an incommensurability in this direction as well. 
Remarkably, the incommensurability A IC , i. e. the displace- 
ment of the peaks from Q AFM , hardly depends on energy 
(see also Fig. 10), in particular there is no indication of a 
sharp, commensurate resonance peak at 70 K. Altogether, 



this indicates that the 70 K data cannot be reproduced by 
a simple temperature-induced broadening of the 5 K data, 
which in turn implies that the dispersion topology is differ- 
ent below and above T c . This is a central point which will be 
further elaborated and confirmed by our discussion of the 
overall dispersion in III D and the fit to the analytical model 
in III E where resolution effects are taken into account. 

Second, the wavevector-resolved representation of Fig. 5 
allows us to conclude that in addition to an intensity redis- 
tribution between different energies, already apparent in 
Fig. 4, there is also an in-plane redistribution at fixed energy 
(Fig. 5g-k; see also Fig. 8e-f). 

Finally, more detailed measurements of the temperature 
evolution of scans through the downward-dispersing branch 
(Fig. 6) show the same abrupt increase in intensity as the res- 
onance peak (panel a and b). 18 ' 80 In addition, panels c) and 
d) show that the onset of superconductivity also abruptly 
changes the geometry, which then remains unaltered upon 
further cooling throughout the superconducting state. 



C. High-energy excitations 

The excitations at high energies differ in several important 
aspects from those at low energies. First, as already men- 
tioned in III A, the in-plane geometry changes from a two- 
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Fig. 5 (color online). Comparison of constant-energy scans at 70 K 
(red circles) and at 5 K (black squares) . The scans along a* are of 
the type (H, -1.5, -1.7), those along b* of the type (1.5, K, 1.7). 
Both datasets were collected under the same experimental condi- 
tions and normalized and BG corrected in the same way. The lines 
show Gaussian fits. The excitation energies denoted in the panels 
are in meV 
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Fig. 6 (color online), a), b) Temperature evolution of the magnetic 
intensity at two points of the downward dispersing mode (co = 
30meV), sketched in reference to the overall hourglass dispersion 
in the SC state. The values were obtained by subtracting a Q-linear 
background, determined by measuring points on both sides of the 
peak, c), d) Q-scan profiles at almost the same fixed energy as in 
a) and b) (co = 31 meV), measured at 5, 50 and 70 K, respectively. 



fold symmetric azimuthal distribution around Q AFM with 
pronounced peaks in the a* -direction below co r to a nearly 
four-fold symmetric distribution with pronounced peaks in 
both diagonal directions above co r (Fig. 3). Whereas the 
detailed model fits of the three-axes data in III E suggest a 
slightly steeper dispersion along b* than along a*, the true 
intensity maxima are between these two directions and do 
not deviate from {100} within the error bars (Figs. 3 and 
13). 

Second, except for a limited energy interval around 
45 meV, where the resonance peak gathers some of its spec- 
tral weight, the onset of superconductivity does not affect 
the spin excitations 17 and in the x"{Q, co) model we will 
use the same function at both 5 K and 70 K. 

Third, whereas we do not detect any intensity at room 
temperature at and below co r , above co r a progressively 
larger fraction of the intensity detected at 5 K is already 
present at 290 K. This is illustrated in Fig. 7, where we 
show scan profiles up to 79 meV From Fig. 4 it is clear that 
this fraction increases mainly because the low-temperature 
intensity decreases, while the high-temperature intensity 
does not change much with energy up to 60 meV 

Finally, we observe that at room temperature the signal 
forms a broad maximum at Qaf M which splits into separate 
peaks away from Q AFM at lower temperatures, but above T c . 
We stress again that the evolving geometry remains four- 
fold symmetric, in contrast to what we see at low energies. 
Overall, the change of the excitations with temperature 
is weaker than at low energies, in agreement with other 
compounds 24 ' 31 ' 108 ' 109 and consistent with the idea of the 
universality of the upper branch. 

A suggestive way to summarize the impact of supercon- 
ductivity on the whole energy evolution is to define four 
energy intervals: The SC-spin-gap interval below ~ 20meY 



the resonance region between ~ 20 and ~ 40meY the 
small "gap" interval around 45meV and the part above 
~ 50 — 55 meV with no superconductivity-induced changes. 



D. Dispersion relations 

Most of the data we show in this paper consist of or is 
based on constant-energy cuts through momentum space. 
There are several ways of representation, each of which 
stresses different aspects of the data. In the previous subsec- 
tions we showed mainly single one-dimensional scan profiles 
or two-dimensional cuts which allow one to see the mea- 
sured data in great detail, but obscure general trends occur- 
ring with changing energy: Overall intensity variations can 
be better visualized in a compact way in three-dimensional 
surface maps of #"(Q, co), like the one shown in Fig. 8. In 
panels a) and b), the large intensity of the commensurate 
resonance peak at co r clearly stands out. However, the evo- 
lution of the Q-pronle away from co r is still not reproduced 
very well. To obtain an even better representation of the 
dispersion, we can take advantage of a redundancy in the 
surface maps: The intensity is coded in both the color and 
the elevation. In Fig. 9 we show the same data as in Fig. 8a)- 
d), but whereas the elevation still represents the intensity, 
the color scale is now normalized to the peak intensity of 
each individual constant-energy scan taken along the black 
lines of constant energy. This relative color scale empha- 
sizes the dispersive character at 5 K in both the a*- and 
b* -direction, with branches merging at the commensurate 
resonance peak 

From Figs. 8c), d) and 9c), d) we see that the hourglass 
shape with commensurate resonance constriction is lost in 
the normal state at 70 K and the intense peak at co r disap- 
pears. Instead, the quasi-one-dimensional signal consisting 
of two Gaussians along a* (see III B) exhibits a nearly verti- 
cal dispersion, forming an overall "Y"-shaped dispersion. 
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Fig. 7 (color online) . Energy evolution of the upward dispersing 
branch in the SC (5 K) and the normal state (290 K). The excitation 
energies denoted in the panels are in meV Data were collected 

in the constant- /c f mode (k { = 4.1 A or 4.5 A ). Due to restric- 
tions arising from the kinematics of the neutron scattering process, 
scans at higher energies were taken in unusually high Brillouin 
zones. They are of the type (H, -2.5, -1.7) and (2.5, K, 1.7) at 
70 meV and (H,H, 1.7) at 79 meV. A Q-linear background has 
been subtracted. The lines result from Gaussian fits to the data. 
The intensities are not comparable among different energies. 
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1(5 K) - 
1(70 K) 




Fig. 8 (color online). Color representation of the mag- 
netic intensity between 26 and 60 meV Panels a), b) 
show the intensity at 5 K, panels c), d) the intensity at 
70 K> T c and panels e),f) the difference between 5 
and 70 K. The upper and lower rows are constructed 
from three-axes scans of the type (H, — 1.5, —1.7) and 
(1.5, iC, 1.7), respectively. The BG was corrected by 
subtracting the intensity at 250 K for co < 38 meV; in 
addition, the data were corrected for a Q-linear BG 
at all energies. The three axes represent transferred 
momentum H or K (r.l.u.), energy (meV) and mag- 
netic intensity (a.u.). Below 38meV, /c f was fixed to 

2.66 A" 1 , above 38meV to 4.5 A~\ Scans at 38meV 
were used to bring both energy ranges to the same 
scale. Crossings of black lines represent measured 
data points. White lines connect the fitted peak posi- 
tions and dotted lines represent upper bounds on the 
incommensurability, as discussed in the text. 



To further elaborate our observations regarding the chang- 
ing signal topology, we show in Fig. 10 the dispersion re- 
lations at 5 and 70 K as obtained from Gaussian fits to the 
scan profiles along a* and b*. Clearly, the signal collapses 
to a single peak at co r in the superconducting state, and the 
plot confirms the finding that the signal exhibits a nearly 
vertical dispersion at 70 K. The evaluation in III E 2 which 
considers the instrumental resolution shows that the data 
are in excellent agreement with a model consisting of two 
incommensurate Gaussian peaks with energy-independent 
A IC below co T and an outward-dispersing damped harmonic 



oscillator mode above co r . However, we cannot exclude that 
at the lowest energies A IC increases slightly, see also Fig. 10. 

In Fig. 10, the right border of the hatched region defines 
an empirical lower bound for the resolvable incommensura- 
bility A IC . While some fits converge towards a A IC within 
the hatched region and thus seem to suggest a persisting 
incommensurability, they provide no significant improve- 
ment over fits with a single peak. This is only partly owed to 
resolution effects 111 - another important reason is that 2A IC 
becomes of the order of the intrinsic width of the peaks. 
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Fig. 9 (color online). Color representation of the magnetic intensity 
between 26 and 60 meX based on the same data as in Fig. 8a) -d), 
but with a different color scale normalized to the peak intensity of 
each individual constant energy scan taken along the black lines. 
Further details are given in the caption of Fig. 8. The panels of this 
figure are reproduced from Ref. 18. 
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Fig. 10 (color online). Dispersion relations along a* and b* at 5 
and 70 K. Here, A IC is obtained from half of the distance between 
the peak positions of two IC peaks. The solid lines are guides to 
the eye. The dashed line represents the dispersion of AF waves 
in insulating YBa 2 Cu 3 6 15 . 110 Incommensurabilities within the 
hatched area come from peaks which lie so close together, that 
they cannot be discriminated from a single Gaussian peak. The 
error in the energy direction is determined by the energy resolution 
and thus depends on k { . Below co r = 38 meV it is of the order of 
2 to 3 mey and above of the order of 6 to 8 meV The error in 
the Q-direction is mainly determined by the statistics of the scan 
(errors of the intensities at the distinct points of the scan) and by 
the question, how well both IC peaks can be separated. From the 
fits we obtain typically < 0.01 r.l.u. in the SC state for large A IC 
and up to ±0.02 r.l.u. at 70 K state close to the hatched region. 
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Fig. 1 1 (color online) . Evolution of constant-energy profiles (raw- 
data) from T = 70 K > T c to T = 5 K in two characteristic energy 
regions, a) Scans along a* at co = co r . b) Scans along b* at 
= 33.5 meV< co r . The data were taken from Fig. 2e) and h). 
At 37meY the data at both temperatures were scaled to the same 
amplitude to allow a better comparison of their Q-widths. Panels c) 
and d) show schematic representations of the dispersion relations 
at 70 K (red) and 5K (dotted grey) in the a*- and b* -direction, 
respectively. Arrows indicate the most salient changes in the SC 
state. 



Hence, at 70 K the scan profiles along b* are indistinguish- 
able from a single Gaussian around and below co r , and we 
can only define an upper limit for a potential incommensu- 
rability, as indicated by the dotted white lines in 8d) and 
9d). Our analysis in III E 2 and the fits in Fig. 17 support the 
interpretation that the signal along b* arises from a single 
peak. 

In Fig. 11c) and d) we illustrate how the resonance forms 
out of the normal-state "Y"-shaped dispersion: Along a*, 
the main effect is the collapse of the incommensurate 70 K 
signal to a single peak, the "neck" of the hourglass dispersion 
(Fig. 11a), accompanied by a moderate signal sharpening 
and increase of A IC below co r . Along b*, on the other hand, 
the main effect is the splitting of the commensurate scan 
profiles into incommensurate peaks below co T , Fig. lib), 
while there is no further sharpening at co r . 

As shown in III A and III C, there is hardly any change 
across T c in the high-energy branch. We thus next exam- 
ine the difference J diff = 7(5 K) - 7(70 K), where 7(5 K) and 
7(70 K) correspond to #"(Q, co) at 5 and 70 K, respectively, 
Fig. 8e) and f). As might be expected, 7 diff is flat at high 
energies and negative around 50 meY see also Fig. 4. At 
low energies, a downward dispersing branch is observed in 
both directions. Due to the subtraction of the commensu- 
rate normal-state signal along b*, the incommensurability 
appears increased in Fig. 8f) . 7 diff thus bears some resem- 
blance to the signal observed in nearly optimally doped 
YBa 2 Cu 3 6 85 and other samples with higher doping, 20 ' 26 
where the lower branch is much more prominent than the 
upper. Regarding the intensity difference is thus very in- 
structive, since it suggests that only the lower branch should 
be associated with the resonance mode. 



E. Analytical formulae for x' '(Q> co) 

Here we will first describe the fitting model for the super- 
conducting state at 5 K (III El) and the normal state at 70 K 
(III E 2) in detail. In III E 3 we compare the models (convo- 
luted with the resolution function) with the corresponding 
experimental data and discuss the limitations of our fit. 



1. Superconducting state (T = SK) 

The best model we found to fit the data at 5 K with a 
reasonable number of parameters is based on a modified 
damped harmonic oscillator. We use the model on a purely 
empirical basis and will not try to justify it physically. 

As suggested by the phenomenology of the data, we sub- 
divide x"(.Q> c°) into an upper (u) and a lower (1) branch, 
both described by a modified DHO function. The upper 
branch describes mainly the the energy range above co r and 
the lower branch mainly the range below co r . However, it 
should be kept in mind that the functions do not drop to 
zero at co r and the intensity of the two branches overlaps. 
It turns out that a smooth cutoff, not contained in the DHO 
model, is necessary to handle these "cross contaminations" 
to the respective other branch. 

The total magnetic susceptibility then reads: 

X "(Q, co) = N [^"(Q, co) + N 0u X ^Q, co)] . (5) 

Here N is an overall normalization factor which is cho- 
sen such that the Q-integrated susceptibility #"(co) has 
a maximal value of 16jtx 2 / (eVf.u.) at the resonance en- 
ergy (Fig. 12), which is known from previous work. 17 ' 80 
N 0u = 1.35 is a scaling factor between the upper and lower 
branch. We first discuss the lower branch (we remind that 
for the sake of a compact notation we have set h = 1 and 
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Fig. 12 (color online). Imaginary part of the local spin suscep- 
tibility x'Xco) for the model functions described in III El (5K) 
and HIE 2 (70 K), compared with the experimental values from 
Fig. 4. The additional points are from Ref. 80 and were obtained 
from measurements on a sample with T c = 67 K, similar to our 
T r = 61 K. 
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Fig. 13 (color online). Constant-energy maps through the model 
for z"(Q> i n the superconducting state at 5 K described in III E 1. 
a) 31 mey b) #"(Q, co) at 31 meV as it would appear in a twinned 
sample. The map was obtained by superposing the map it a) 
with its transpose, c) 60meV The corners of the black square 
define the positions of maximal intensity of the upper branch in 
the {110} directions. The shown intensity was not convoluted with 
the instrumental resolution function, is shown in arbitrary units 
and is not comparable between different panels. 



in addition there is no anisotropy in the damping or the 
amplitude (Fig. 13). The pref actor N u (co,Q) is somewhat 
more complicated than A^(co), since it accounts for three 
effects: the cutoff towards low energies, N^^co), the 
superconductivity-induced suppression above co r , iV dip (co) 
(Fig. 12) and the azimuthal intensity distribution with peaks 
along the (110) directions, N rot (Q) (Fig. 7), 

N u (co, Q) = iV cutoff (co)iV dip (co)iV rot (Q). (12) 

The different terms are 

1 

NnitnffC&O = 1 T, ~ ~, (13) 

cutofa J l + exp[(co-36)/1.5] 

N dip (co) = 1 - 0.2 exp [-{co - ft dip ) 2 /(2a 2 ip )] , (14) 



quote co in meV) 

xl'(Q,o>) = N t (a>) 



2C0redr Q 



fared- +( 2 ^redr Q y 



(6) 



where co red = 2co r - co and co r = 38.5 meV is the resonance 
energy. 112 The dispersion of the mode is given by 



co Q = co r + s G (H - H ) 2 + s b (K- K )\ 



(7) 



where Q AFM = (H ,K ) is the antiferromagnetic wave vector 

o 2 o 2 

in the plane, with s a = 280meVA and s b = 455meVA . 
The damping y Q has the following angular dependence: 



r Q = r a + Ay 



(K-K f 

Q ~~ QafmI 2 



(8) 



Here a denotes the a* -direction in the Brillouin zone and 
Y a = 1.5 meV and Ay = y b — y a = 6.5 meV 113 The prefactor 
Ni(co) implements the superconductivity-induced spin gap 
below ~ 24 meV and the mentioned cutoff at high energies 
above ~ 42meV: 114 



N,(co) = 



exp [(co - 24)/3] 



1 



1 + exp [(co - 24)/3] l + exp(co-42) 
The upper branch is given by 

2 cor 

Xu(Q,co)=N u (to,Q) 



(co 2 -^ 2 ) +(2cor) 2 



(9) 



(10) 



where in contrast to the lower branch the damping r = 
11 meV is a constant, independent of the azimuthal angle. 
Also, due to the steeper dispersion at higher energies, we 
model ft Q with a quartic power law (p = 4) : 



= co r + [s 2 Jp (H - H ) 2 + S 2 b /P (K - K f] ?/2 . (11) 



Here S a = 4830meVA P and S a = 10065 meVA P . Despite 
the seemingly large difference between S a and S b , the in- 
plane anisotropy is smaller than in the lower branch, and 



iV rot (Q) = 1-0.3 



(H-H f-(K-K f ^ 2 
IQ- QafmI 2 



(15) 



and n dip = 47 me\J cr dip = 2meV The whole function 
x"iQ> °>) needs to be cut at the zone center and contin- 
ued periodically to account for the lattice symmetry. 



2. Normal state at 70 K 

At 70 K we apply the same subdivision into a lower and 
upper branch like in eq. (5) . The prefactor N has the same 
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Fig. 14 (color online). Color representation of the magnetic in- 
tensity between 26 and 60 meV at T = 5 K. a), b) Measured data, 
identical with Fig. 9a), b). c), d) Fit to the damped harmonic 
oscillator model described in eqs. (5)-(15) of HIE 1. A detailed 
comparison for the individual scans is shown in Fig. 17. 
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Fig. 15 (color online). Color representa- 
tion of the magnetic intensity between 
26 and 60 meV at T = 70 K. a), b) Mea- 
sured data, identical with Fig. 9c), d). 
c), d) Fit to the model described in 
III E 2, which best reproduces the data. 
A detailed comparison for the individual 
scans is shown in Fig. 17. e), f) Repre- 
sentation of a broadened damped har- 
monic oscillator model (increased damp- 
ing constants) which fails to reproduce 
salient features of the data, like the ab- 
sence of an hourglass constriction and 
the nearly vertical dispersion. 



numerical value like at 5 K. The upper branch x'^iQ, co) is 
described by the same model (10) as at 5 K, with the same 
dispersion (11) and damping r = 11 meV Also the pref actor 
N rot (co) is the same as in (15), while iV dip being 1 and iV cutoff 
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Fig. 16 (color online). Color representation of the subtracted mag- 
netic intensity 7(5 K) - 7(70 K) between 26 and 60meV a), b) 
Measured data, identical with Fig. 8e), f). c), d) Difference be- 
tween the DHO model describing the data at 5 K, III E 1, and the 
model describing the data at 70 K, IIIE2. 



being 



cutoff 1 



(co)- 



1 



1 + exp [(co - 37)/2] 



(16) 



are different from (14) and (13), respectively. 

In contrast, no satisfactory fit can be obtained for the low- 
energy branch at 70 K with the DHO model used at 5 K. The 
best one can do within the framework of the DHO model is 
to increase the damping parameters y a and Yb significantly 
- however, the hourglass character of the dispersion is not 
lost by this (Fig. 15e, f). A much better fit is obtained 
if the lower branch is modeled by two anisotropic two- 
dimensional Gaussian distributions, displaced from Q AFM 
along the a* -direction: 



i=l,2 



-4 In 2 



f H ~ Hi Y | ( K ~ Ki 



(17) 



Within the energy range of ~ 25 — 38 meV the experimen- 
tally observed incommensurability does not change within 
the error, hence we can assume Q 1 = (H 1 ,K 1 ) = (0.575, 0.5) 
and Q 2 = (H 2 ,K 2 ) = (0.425,0.5) to be energy-independent. 
The Gaussian widths are a a = 0.125 and a b = 0.17, both 
in r. 1. u. The normalization factor is given by 



N^co) = No^N^pfcojN^td^fcojN^fco), 
where N oi = 0.0181. 



NcutupO) = 



1 



1 + exp [(co - 37)/2] 

provides a cutoff towards higher energies, while 

_ 0.5 + exp[(co-29.5)/2] 
Vcutdown(^) - 1 + exp [(co-29.5)/2] 



(18) 



(19) 



(20) 
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Fig. 17 (color online). Comparison between the measured data (as shown in Figs. 5-8, full squares) and the corresponding model (open 
circles): The DHO model, for T = 5K (HIE 1, left two columns); and the model for T = 70 K (HIE 2, right two columns). 
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Fig. 18 (color online). Continued from Fig. 17: Comparison be- 
tween measured data (full squares) and the DHO model, III E 1, at 
T = 5 K, (open circles). 



provides a partial cutoff towards lower energies. Finally, 

N ^=\co/27 forco<27meV (21) 

provides a simple linear decrease below 27meV consistent 
with previous work. 106 



3. Summary of the modeb 

In Fig. 12 we show the momentum-integrated spin 
susceptibility ^"(Q, co) for both models (5K and 70 K) 
in absolute units. The sum-rule integral S = 
(TigVg)" 1 J Q dtodQx"(.Q,co) is identical for both mod- 
els within the error bars, as required. The fact that the 
experimental Q-integrated spin susceptibility from Fig. 4 
agrees well with the model curves lends confidence to our 
results. The comparison with experimental values previ- 
ously obtained 80 in a sample with a T c of 67 K (Fig. 12) 
indicates that our extrapolation towards high energies is 
correct up to at least 100 meY 115 and with a somewhat in- 
creased error tolerance probably far above these energies. 
This is further supported by a comparison with the evolu- 
tion of the high-energy data of Hayden et al. 17 The in-plane 
maps presented there suggest that the in-plane intensity 
distribution of our model, in particular the axial intensity 
distribution with maxima in the {110} directions, agrees 
with experiment up to the same energy. This is also apparent 
from a comparison with Fig. 13 where we show constant- 
energy maps of #"(Q, co) at energies below and above co T . 
When comparing it with Fig. 3 and other measured data, it 
should be kept in mind that Fig. 13 shows intensity that was 
not convoluted with the instrumental resolution function 
or averaged over a certain energy range: In the measured 
data, the intensity distribution always appears spread over 
a larger region. We also mention that due to the decreas- 
ing intensity and spurious contaminations of the signal it is 
difficult to obtain detailed data far below 20 meV We can 



clearly establish the opening of a superconducting gap and 
find indication for a linear decrease of the intensity. Also the 
strongly anisotropic incommensurability we observe at the 
lowest covered energies implies a natural extension towards 
co = 0. All this is in agreement with Refs. 80 and 106 and 
suggests the validity of our models for energies approach- 
ing co = 0. Finally, we also note that our model (Fig. 13b, 
artificially "twinned") is in remarkable agreement with the 
"castle"-like structure reported by Mook et al. (Ref. 97) for 
a twinned YBa 2 Cu 3 6 6 sample. 

We will now evaluate in more detail how well our models 
reproduce the experiments. To do this, we compare the 
model functions, convoluted with the instrumental resolu- 
tion function, with the experimental data using two different 
representations: In Figs. 14-16 we use the color representa- 
tion introduced in III D (Figs. 8 and 9), which permits the 
assessment of the salient features like dispersion, hourglass 
constriction and resonance peak at a glance. Figs. 17 and 
18, on the other hand, allow us a detailed comparison of 
the model with the individual constant-energy scan profiles. 

Beginning in the superconducting state at 5 K, the two rep- 
resentations show the excellent fit of our DHO model, both 
qualitatively and quantitatively. In particular, the hourglass 
dispersion with the typical intensity increase and Q-width 
constriction at co r is reproduced very well. The comparison 
in Fig. 17 shows that this is not only true for the shape but 
also for the intensity. In Fig. 18 we compare our model 
with a scan profile along a* which does not cross Q AFM but 
Q = (0.5, 1.44, 1.7). The good agreement confirms that our 
model also captures the essential features of the detailed 
in-plane distribution of the lower branch, as was already 
indicated by a qualitative comparison of Fig. 3a,b) and 
Fig. 13a). 

Proceeding to 70 K, we first note that the data cannot be 
reproduced by simply broadening the DHO model function 
used at 5 K. This is illustrated by a comparison of the data 
(Fig. 15a,b) with the "best" fitting broadened DHO model 
(Fig. 15e,f). In particular, the DHO model still exhibits a 
commensurate resonance constriction, and the the lower 
and upper branches disperse towards the constriction in 
a way comparable to the situation at 5 K. The only well- 
reproduced part of the spectrum is the high-energy part, 
which does not change significantly compared to 5 K anyway. 

A much better agreement with the data is obtained using 
the two-Gaussian model described in III E 2 for the lower 
branch: The steep dispersion and the absence of an hour- 
glass constriction are reproduced very well. Also the quasi- 
one-dimensional character at co < co r is captured. The de- 
tailed comparison with the individual scan profiles (Fig. 17) 
again confirms the good agreement. Both representations 
show that the signal along b* remains commensurate down 
to lowest covered energies. 

Next, we show in Fig. 16 that our models also reproduce 
the difference between the spin susceptibilities at 5 and 
70 K, where we observe just a downward dispersing mode. 

Finally, we revisit a point raised in III D, namely the dis- 
persion relations of the upward- and downward-dispersing 
modes. In Fig. 19 we compare the dispersion obtained from 
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Fig. 19 (color online). Comparison of the dispersion relations as 
obtained from Gaussian fits to the data (symbols, see Fig. 10) with 
those of the model derived in III E 1 for T = 5 K (solid lines). Black 
symbols and lines are for the a*-direction in a) and the b*-direction 
in b) . Stars and dash-dotted lines are for the diagonal direction, 
where A IC was obtained from the TOF data for the points at 63 
and 75 meV (Fig. 3), from the 3-axis data at 79 meV (Fig. 7) and 
from TOF measurements shown in Ref. 17 for 85 meV 



Gaussian fits to the scan profiles with the dispersion of the 
model functions. An inspection of Fig. 13 explains why the 
Gaussian fits always tend to yield a smaller incommensura- 
bility A IC than the model #"(Q, co): The in-plane azimuthal 
distribution exhibits appreciable intensity not only in the 
principal crystallographic directions but also between them. 
Due to the large resolution volume, in the course of a scan 
along a* or b* (or along {110}, for that matter) intensity 
is still collected far beyond the maximum of ^"(QjCo) in 
this particular direction. This shifts the intensity maximum 
of the measured scan closer to Q AFM , in particular at high 
energies where the resolution volume is especially large 
and the azimuthal intensity maxima are along {110}, af- 
fecting the scans along a* and b*. This demonstrates that 
it is safer to consider the full in-plane distribution of the 
model x"(Q, co), extracted under inclusion of the resolution 
function, as we did. Just showing the intensity maxima in 
particular directions, no matter whether extracted from sim- 
ple fits to the data or from the model #"(Q, co), reproduces 
qualitative features but only conveys limited quantitative 
information. 



IV DISCUSSION 

In the previous section III we have shown that the salient 
features of the spin excitations in YBa 2 Cu 3 6 6 do not arise 
from resolution effects. Our analysis in III E has led to quan- 
titative models for #"(Q, co) in both the superconducting 
and the normal states which capture these features, namely 
the continuous low-energy excitation branch at 5 K and the 
abrupt topology change towards a "Y"-shaped dispersion 
with no pronounced resonance feature at 70 K; the change 
of the spectral distribution at low energies from weakly to 
strongly anisotropic upon heating above 70 K; and the four- 
fold symmetry of the upper branch and its insensitivity to 



temperature variations across T c . 

What are the implications of these results for the scenarios 
invoked to explain the peculiar spin excitations and the 
delicate interplay between different competing phases in 
underdoped cuprates? 

When the hourglass dispersion was reported for the 
La 2 _ x (Sr,Ba) x Cu0 4 family, 22 the notion of bond-centered 
stripes in the form of weakly-coupled two-leg spin ladders 
was put forward as an explanation, 22 and the concomitant 
observation 17 of a similar dispersion geometry in the super- 
conducting state of YBa 2 Cu 3 6+JC seemed to establish this 
as a universal scenario for cuprates. This idea has triggered 
many subsequent theoretical papers 38-41 but it has now be- 
come clear that it is not applicable to YBa 2 Cu 3 6+x without 
further elaboration. 18 The first difficulty is the observed 
change from two-fold symmetry at low energies to four-fold 
symmetry at high energies (Figs. 3, 7 and 13). In the model 
of two-leg spin ladders, a constant-energy cut through the 
upper branch should show two parallel streaks running 
across the ladder, since the gapped excitations disperse only 
along the one-dimensional ladder. Introducing weak corre- 
lations between the ladders results in two incommensurate 
spin-wave-like cones at low energies, displaced from Q AFM 
perpendicularly to the ladder direction. The resulting spec- 
trum is two-fold symmetric at any energy, in contrast to our 
data from twin-free samples, which show incommensurate 
peaks along both a* and b* in high-energy scans through 
Q AFM , (Fig. 7). One might object that orthorhombicity does 
not provide a sufficiently strong potential to pin ladders 
either along a* or along b*, and that domains of both orien- 
tations are present. This scenario is implicitly included in 
the calculations of Refs. 38-41, where the simulated data 
are artificially twinned. However our observation of a strong 
in-plane anisotropy at low energies in the normal-state data 
contradicts this assumption. 

Different approaches have been used to address and to 
some degree to resolve these problems. Magnetic fluctu- 
ations of unidirectional, metallic stripes numerically com- 
puted in the framework of a Hubbard model exhibit two-fold 
symmetry at low energies, and high-energy branches dis- 
persing along both orthogonal principal crystallographic 
directions without the necessity to invoke twinning. 42 How- 
ever, the model predicts a gap between lower and upper 
branches in the a* -direction and a non-monotonic evolution 
in the b* -direction, which we do not observe experimentally. 
It has also been shown that the "Y"-shaped dispersion can 
be understood as arising from fluctuating short-range stripe 
segments, and that an unequal population of the two perpen- 
dicular stripe orientations results in increasing anisotropy at 
low energies, in agreement with our observations. 52 How- 
ever, this model is purely phenomenological and contains 
many fitting parameters that cannot be determined from 
first principles. 

A second shortcoming of these early models of the hour- 
glass dispersion is their failure to address the influence of 
superconductivity on the spin fluctuation spectrum. The 
neutron data we have presented here clearly document a 
transition from a "Y"-shaped to an hourglass dispersion at 
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T <T C . We note that this change of the dispersion topology, 
and in particular the fact that the spectrum at 70 K cannot 
be obtained from the one at 5 K by increasing the mode 
damping (HIE 2 and Fig. 15), is difficult to reconcile with 
the notion of a collective mode above T c , 116 which is merely 
concealed due to scattering by continuum excitations. 

Our data rather indicate a "resonant" collective mode 
that only appears below T c . The occurrence of a resonant 
mode in the d-wave superconducting state has been the- 
oretically treated by many authors/ 1-78 ' 117 as reviewed in 
Ref. 53. In a simple Fermi-liquid-based scenario, the reso- 
nance is described as a spin exciton, i. e. a spin-triplet bound 
state stabilized by electron-electron interactions below the 
electron-hole continuum setting in at twice the supercon- 
ducting gap. The observed intensity redistribution below 
T c , in which the resonant mode draws spectral weight from 
both lower and higher energies (Figs. 8e),f) and 16), as 
well as the formation of a downward dispersing branch 
with a shape roughly following the d-wave superconducting 
gap support this kind of scenario. 118 In the framework the 
weak-coupling Fermi-liquid scenarios proposed so far, 119-121 
hybridization of electronic states based on the CuO chains 
and Cu0 2 layers in the YBa 2 Cu 3 6+x crystal structure with 
electronic states in the Cu0 2 layers might provide an expla- 
nation of the modest in-plane anisotropy in the spectra in 
the superconducting state, but appears too weak to account 
for the strong anisotropy of the spectra above T c (Fig. 5) . We 
also note that the results of most weak-coupling calculations 
have thus far been presented in arbitrary intensity units. In 
view of the recently reported failure of such calculations 
to reproduce the absolute intensity of the observed spin 
excitations in an electron-doped cuprate superconductor, 122 
the comparisons between experimental and numerical data 
on an absolute intensity scale constitutes an important chal- 
lenge for the future. The results presented here provide a 
good basis for studies designed to meet this challenge. 

Recent theoretical work indicates that it may be possible 
to describe the spin excitation spectra above and below T c 
in the framework of a single scenario, according to which 
superconductivity competes with a Pomeranchuk instability 
of the Fermi surface. It is noteworthy that from a symmetry 
point of view, the Pomeranchuk model and the fluctuating- 
stripe models are equivalent: In both cases, the C 4 symmetry 
is spontaneously broken. An analysis of #"(Q, co) in the 
slave-boson mean-field scheme of Ref. 36 reproduces the 
transition from a "Y"-shaped dispersion to an hourglass 
dispersion upon entering the SC state and (qualitatively) 
also the crossover from two-fold to four-fold symmetry with 
increasing energy. Future work will have to show if this 
framework can be extended to even lower doping levels, 
where the ground state exhibits incommensurate magnetic 
order and nearly gapless magnetic excitations, 31 ' 34 and is 



believed to become insulating at high magnetic fields. 123 
In this regime, the localized-electron scenarios that have 
been discussed in the context of the spin excitations of 
La 2 _ JC (Sr,Ba)j C Cu0 4 may be a more promising starting point. 

For a complete understanding, some questions remain to 
be addressed in future investigations: First, it is interesting 
to determine the dispersion topology in YBa 2 Cu 3 6 45 and to 
compare it to the typical "Y"-shaped dispersion we see here. 
Second, the onset temperature of the incommensurability 
in YBa 2 Cu 3 6 6 must be determined. This is particularly 
important in light of the fact that, based on the onset tem- 
perature of the in-plane anisotropy of the Nernst effect, T v 
the authors of Ref. 63 relate electronic liquid-crystalline 
order to the pseudo-gap phenomenon. Interestingly, the 
onset temperature of the novel magnetic order at Q = 0, 
which was also related to the pseudo-gap state, is ~ 220 K 
in our YBa 2 Cu 3 6 6 sample, 27 and thus close to T v . 

Finally, we expect that the full Q- and co -resolved data set 
of magnetic excitations and the analytical description we 
have presented here will greatly facilitate quantitative com- 
parison with spectroscopic data obtained by complementary 
experimental methods. Previous attempts to relate anoma- 
lous features in charge -fluctuation spectra to the spin exci- 
tations have often narrowly focused on particular regions of 
momentum space. For instance, because of the strong tem- 
perature dependence of the spin fluctuation spectrum near 
Q AFM , the insensitivity of the nodal kink in ARPES spectra 
to changing the temperature across T c had been taken as 
an indication that this feature originates from coupling to 
phonons. We have shown, however, that coupling to the 
weakly temperature dependent high-energy spin-excitation 
branch contributes substantially to the nodal kink observed 
in ARPES spectra. 33 For the same reason, the absence of 
renormalization effects attributable to the magnetic reso- 
nant mode in optical spectra of overdoped cuprates 124 does 
not invalidate theories based on magnetically mediated pair- 
ing mechanisms. A quantitative description of optical and 
STS spectra in terms of spin-fermion models requires a full 
energy- and momentum-integration over all spin-fluctuation 
mediated scattering channels. The analytical model func- 
tions presented here for YBa 2 Cu 3 6 6 provide a good start- 
ing point for such endeavors. As we have shown in the case 
of ARPES, 33 this approach also allows quantitative experi- 
mental tests of spin-fluctuation mediated pairing models of 
high-temperature superconductivity. 
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